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Abstract. The defect of a function f : M R is defined as the 
difference between the measure of the positive and negative regions. In 
this paper, we begin the analysis of the distribution of defect of random 
Gaussian spherical harmonics. By an easy argument, the defect is non- 
trivial only for even degree and the expected value always vanishes. Our 
principal result is obtaining the asymptotic shape of the defect variance, 
in the high frequency limit. As other geometric functionals of random 
eigenfunctions, the defect may be used as a tool to probe the statistical 
properties of spherical random fields, a topic of great interest for modern 
Cosmological data analysis. 



1. Introduction 

In recent years, a lot of interest has been drawn by the analysis of the 
geometric features of random eigenfunctions for the spherical Laplacian. 
More precisely, let us denote as usual by {^m(-)}m=-« ; , for / = 1,2,... 
the set of (real-valued) spherical harmonics, i.e. the array of real-valued 
functions satisfying 

As2Yim = + l)Yim, 

where A_52 denotes the spherical Laplacian (see for instance [T7], [Ej). The 
random model we shall focus on is 

1 ' 

(1) fi{x) = /^TTT ^ 

m=—l 

where the coefficients {aim} are independent standard Gaussian with zero- 
mean and unit variance. 

The random field /; is isotropic, meaning that for any S N and xi, . . . , G 
5^, the distribution of the random vector {fi{xi), . . . , fi{xi)} is invariant un- 
der rotations, i.e. equals to the distribution of {fi{g • xi), . . . , fi{g ■ xi)} for 
any rotation g E SO (3) on the sphere. Also, /; is centred Gaussian, with 
covariance function 

(2) ri{x,y) := E[Mx) ■ fi{y)] = Pi {cos {d{x,y))), 

where Pi are the usual Legendre polynomials defined by Rodrigues' formula 
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and d{x, y) is the (spherical) geodesic distance between x and y. As weh- 
known, Legendre polynomials are orthogonal w.r.t. the constant weight 
uj{t) = 1 on [—1, 1], see Section [5] for more discussion and details. 

1.1. Background. Random eigenfunctions for the spherical Laplacian nat- 
urally emerge in a number of different physical contexts. A particularly ac- 
tive area is related to the analysis of isotropic spherical random fields on the 
sphere, as motivated for instance by the analysis of Cosmic Microwave Back- 
ground radiation (CMB), see for instance O [6] and the references therein. 
Under these circumstances, eigenfunctions like fi represent the (normalized) 
Fourier components of the field, i.e. the following orthogonal expansion 
holds, in the sense: 

oo 

f{x) = ^cifi{x) , 

1=1 

where q > is a deterministic sequence (depending on the angular power 
spectrum) which encodes the full correlation structure of /. 

In the CMB literature, f{x), and hence the components fi, are actually 
observed by highly sophisticated satellite experiments such as WMAP by 
NASA and Planck by ESA. It is then common practice to analyze geomet- 
ric functionals of the observed CMB radiation to constraint the statistical 
properties of the underlying fields, e.g. to test for isotropy and/or Gaus- 
sianity. For instance, the three Minkowski functionals, providing the area, 
the boundary length and the topological genus of excursion sets over a given 
level {A^ := {x £ S"^ : f{x) > z}) have been applied on CMB data by a huge 
number of authors, including [8l[9l[13]. Many efforts have been spent to ana- 
lyze unexpected features; for instance, the so-called Cold Spot (see Oil] and 
the references therein), whose statistical significance can also be evaluated 
by means of local curvature properties of isotropic Gaussian fields ( [7J ) . 

In this work, we shall focus on one of the most important geometric 
functionals, namely the defect. The defect (or "signed area", see [2]) of a 
function -0 : — t- M is defined as 

(3) P(V') := meas (^^"^(0,00)) - meas (V'"^(-oo,0)) = / n{ip{x))dx. 
Here 7i{t) is the Heaviside function 

(4) ^(t) = l[o,oo)(i)-l(-oo,o](*) = 



1 t > 
-1 t < 
t = 



where lA(i) is the usual indicator function of the set A, and dx is the 
Lebesgue measure. The defect is hence the difference between the areas of 
positive and negative inverse image of ip, respectively. 

Of course, the defect Vi = V{fi) of fi is a random variable; its distribu- 
tion is the primary focus of the present paper. Note that for odd /, // is 
always odd, so that in this case the defect vanishes identically (P; = 0), and 
therefore T>i has nontrivial distribution for even I only. A possible alterna- 
tive to avoid trivialities is to restrict ourselves to a subset of the sphere, the 



THE DEFECT VARIANCE OF RANDOM SPHERICAL HARMONICS 



3 



most natural choice being a hemisphere, i.e. we may choose any hemisphere 
£ and define 

Vf := ju{fi{x))dx- 
in this paper we study T>i for even / only. 

1.2. Statement of the main result. Evaluating the expectation of Vi is 
trivial. Indeed, notice that integration over 5^ is exchangeable with expec- 
tation, so that 



(5) Var(PO = 72(1 + 0(1)), 



and K[T-L{fi{x))] = vanishes for every x G 5^, by the symmetry of the 
Gaussian distribution. We just established the following lemma: 

Lemma 1.1. For every I = 1,2, we have 

nvi] = 0. 

The main result of the present paper concerns the asymptotic behaviour 
of the defect variance: 

Theorem 1.2. As I — >• 00 along even integers, the defect variance is as- 
ymptotic to 

C. 
P 

where C > is a positive constant. 

The constant C in (l5|) may be expressed in terms of the infinite (condi- 
tionally convergent) integral 

00 

(6) C = 32tt j ip (arcsin(Jo(V')) - Jo(V')) d^- 



See Section [2] for some details on the integral on the RHS of dH) . We do not 
know whether one can evaluate C explicitly; however, we shall be able to 
show that 

c > 



(see Lemma lS.Sp . 

One should compare the statement of Theorem 11.21 to empirical results of 
the study conducted by Blum, Gnutzmann and Smilansky The authors 
of that work studied the defect (or, as they refer to, the "signed area") of 
random monochromatic waves on various planar domains. For the partic- 
ular case of unit circle, they found that the order of magnitude of defect 
variance is consistent to (0) per unit area, with leading constant evaluated 
numerically as ~ 0.0386. The corresponding defect variance per unit area 
in our situation is 



with leading constant 



C > = = 0.0389... 

7r2 X y/27 
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1.3. Previous work. To put our results in a proper perspective and explain 
the technical difficulties to be handled, we need to briefly recall some results 
from [12]. In that work we studied the asymptotic behaviour of excursion 
sets 

of as / — )• oo. As we will recall below, a rather peculiar phenomenon can 
be shown to hold for z ^ 0; namely, the asymptotic distribution of the area 
of these excursion sets is fully degenerate over z, i.e. it corresponds to a 
Gaussian random variable times a deterministic function of z. This result is 
intuitively due to an asymptotic degeneracy in the excursion set functionals, 
which turns out to be dominated by a single polynomial (quadratic) term. 

This degeneracy, however, does not hold for the special case z = 0. More 
precisely, in [12] we focused on the asymptotic behaviour of the empirical 
measure of random spherical harmonics, defined as 

^l{z) := / li^ c<,)ifl{x))dx = meas{x : fi{x) < z} , 

z gM. A key step in that paper is the asymptotic expansion 

(7) <^i{z) = 47r X cD(z) + V ^hi,g, hi,, := / H,{fi{x))dx, 

where Hq[.) are standard Hermite polynomials, ^{z) = Pr{Z < z} is the 
cumulative distribution function of a standard Gaussian variable, and the 
deterministic functions Jq{z) can be explicitly provided in terms of higher 
order derivatives of ^{z), Jq{z) = {—lY<^^'^\z) (see [l^ for more discussion 
and details). 

It turns out that, as / — t- oo, 

^hi,2 A^(0, 1) , Vlhi.q = Op(l) , for g > 3; 

moreover J2{z) = —z(j){z) clearly does not vanish for all z ^ 0, whence the as- 
ymptotic behaviour of Vl{^i{z) — 4it x ^{z)} is easily seen to be Gaussian, 
uniformly over z. Furthermore, the limiting process is completely degenerate 
with respect to z, a feature to which we shall come back later. 

For 2; 7^ 0, the asymptotic behaviour of the area functional for the excur- 
sion sets is hence fully understood. The previous argument, however, fails 
for 2: = 0, as in this case the leading term is null and each summand in the 
asymptotic expansion ([7]) becomes relevant. Up to a linear transformation, 
this is clearly equivalent to the defect functional, indeed 

Thus, the case z = is the most challenging from the mathematical point 
of view, and, at the same time, the most interesting from the point of view 
of geometric interpretation. 

1.4. Overview of the paper. The plan of the paper is as follows. In 
Section [21 we provide the main ideas behind our principal arguments, to 
help the reader understand the material to follow; in Section [3l we discuss 
the relation of our results to recent works on the distribution of nodal lengths 
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and level curves for random eigenfunctions and related conjectures; Section 
Hlprovides the proof of the main results, whereas Section [5] contains auxiliary 
lemmas, we believe of independent interest, on the asymptotic behaviour of 
moments of Legendre polynomials. 
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2. On the PROOF OF Theorem IL2J 



To establish our results, we shall need a detailed analysis of the odd 
moments of Legendre polynomials 



7r/2 



Pi (cos I 



\2fe+l 



sm 



The rationale for this can be explained as follows. It is relatively easy to 
express (up to a constant) the defect variance as 



7r/2 



(8) 



arcsin(P;(cos 9)) sin 9d9, 



where sin 9d9 is (up to a constant) the uniform measure on the sphere in the 
spherical coordinates (see Lemma l4.ip . It then remains to understand the 
asymptotic behaviour of a task which is put forward in Proposition 14. 2| 
here we provide the main ideas underlying its proof. We know from Hilb's 
asymptotics (Lemma 15. ip . that P;(cos6') has a scaling limit: as / — ?■ oo, for 
any fixed 



we have 



(9) 



Pi I cos 



0,(/ + l/2 



vr 

'2J 



/ + 1/2 



Jo(V')- 



In fact, the latter estimate holds uniformly for ip = o{l), and it may be 
shown that, as a tail of a conditionally convergent integral, the contribution 
of the other regime tp ^ I is negligible (see the proof of Proposition 14. 2p . 

In this sequel we neglect the difference between / and I + ^- It is then 
natural to try to replace P;(cos0) in ([8]) with its scaling limit; a formal 
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substitution yields heuristically 
(10) 

«7r/2 I-k/2 

f« y j arcsin {Pi {cos {ip / 1))) sin (ip / 1) dip ^ J J arcsin(Jo(V')) sin (-(/;//) d-^ 


Itt/2 

~ ^ y arcsin(Jo(V'))V''^V') 



where we replaced sin with j, which is justified for tp = o{l). This is 
inconsistent to the statement of Theorem \1.2[ the integral 

oo 

(11) / arcsin( Jo(V'))V''^V') 



diverges, so some more care is needed to transform (jlOp . Moreover, the last 
integrand in (|lUp is somewhat different from the integrand in the integral 
dU]) defining C (neglecting the constant in front of the integral in ([U])) in that 
in 1^ we subtract Jo(V') from the arcsine. Note that the factor ip in both 
(jlOp and ^ is reminiscent of the uniform measure sin 9d9 on the sphere. 

There is a subtlety that explains the latter discrepancy between the in- 
tegrand in the definition ([HD of C and the integrand in (fTUj) . One way to 
justify the substitution of the scaling limit of Pi{cos9) (i.e. the second step 
of (jlOp ) is expanding the arcsine in the integral of ([8]) into the Taylor series 
around the origin 

oo 

arcsin(t) = afct^^'+\ 

k=0 

for some explicitly given coefficients (see ()19p and the formula immedi- 
ately after). We then need to evaluate all the odd moments of Pi{cos 0) w.r.t. 
to the measure sin 6 dO on [O, It turns out that the moments (appropri- 
ately scaled) are asymptotically equivalent to the corresponding moments of 
the Bessel function on M-|_ w.r.t. the measure ipdip, with the exception for 
the first moment (Lemma l5.2p . Indeed, Pi{cos9) integrates to zero, whereas 

oo 

J Jo{ip)'ipdip diverges; this phenomenon also accounts for the divergence of 



(jlip . To account for this difference, we need to subtract Jo{ip) from the 
arcsine in the integrand (|lUp . thus obtaining a conditionally convergent in- 
tegral. This is indeed the heuristic explanation for the discrepancy between 
(USD and ©. 

To summarize this discussion, we note the following. The variance is 
given, up to an explicit constant, by the integral dH). The constant C in ([6]) 
is reminiscent of // up to the scaling ([9]): the expression Pi{cos 9) is replaced 
by its scaling limit Jo{ip), and the measure sin 9d9 becomes ipdip on after 
scaling. Subtracting Jo(V') from the integrand in Q accounts for the fact 
that for even I, the integral of P/(cos^) against sin^d^ vanishes, whereas 
this is not the case of integral of the scaling limit Jo(V') against ipdip. 
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3. Discussion 

As reported in Section 11.31 the asymptotics for the variance of the excur- 
sion sets is given for z / by 

Var($;(z))~z2^(z)2.— 

(here (j) is the Gaussian probabihty density); for the defect (z = 0), this 
gives only an upper bound o (j). The principal result of the present paper 
(Theorem II. 2|) states that the defect variance ^ is of order of magnitude 
Our explanation for this discrepancy is the disappearance, for z = 0, of 
quadratic term in the Hermite expansion of the function l.{fi < z), see |12| . 

We should compare this situation to the length distribution of the level 
curves. For t G R let 

4 = £*(/,) = length ifrHt)). 

be the (random) length of the level curve /^"^(t), the most important case 
being that of i = 0; the corresponding level curve is called the nodal line. It 
is known [20] that for t S M, the expected length is 



E [Cj] ~ cie-*'/2 . i_ 
The variance is asymptotic to 

Var (4) ~ C2e-*'t^ • / 
for t ^ 0, whereas for the nodal length 

Var(£[') ~C3log/. 

It is therefore natural to conjecture that the discrepancy of the defect and the 
empirical measure are somehow related to the discrepancy in the variance 
of the nodal length. However, no simple explanation, as in the case of the 
defect, is known to explain the discrepancy between the nodal length and 
the non-vanishing level curves length. 

Note that we may relate to the excursion sets by the identity 



(12) $/(z) = J C\dt. 

— oo 

Let us recall that for two random variables X,Y , the correlation is defined 
as 

Cov{X,Y) 



Corr(X, Y) 



V^Var(X)^Var(y) 



It is known [2D] that for every ti,t2 G become asymptotically fully 

dependent for large I, in the sense thalQ 

(13) C0TT{C\\Cf)=l-0i^^{l). 
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In fact, a much stronger statement regarding the rate of convergence was 
proven. It then seems, that the asymptotic degeneracy of ^i{z) is, essen- 
tiahy, an artifact of (jl2p and the asymptotic full dependence of the level 
lengths (fT3ll . 

A possible explanation for the phenomenon ()13p is the following conjecture 
due to Mikhail Sodin. Let x £ and for t G M let Cj.^ = CUfi) (the "local 

length") be the (random) length of the unique component of fi~^{t) that 
contains x inside (or 0, if /; does not cross the level t). 

Conjecture 3.1 (M. Sodin). The local lengths are asymptotically fully de- 
pendent in the sense that for every x G 5^ and ti,t2 S M 

Corrf4^,4^U 1-0,^00(1). 



''l;x ' l;x 

Heuristically, Conjecture 13.11 is stronger than (and thus the asymp- 
totic degeneracy of the empirical measure via p2p ). since can be viewed 
as a summation of £^ over a set of points x E 5^. More rigorously, we 
should assume that fi{x) 7^ t; it is immediate to see that this is satisfied 
almost surely. 



4. Proof of Theorem 11.21 

The first Lemma is probably already known, but we failed to locate a 
direct reference and thus we report it for completeness. The proof is based 
on standard properties of Gaussian random variables. 

Lemma 4.1. For I even we have 

7r/2 

(14) Var(Pz) = 327r J arcsin(Pi (cos (9)) sin (9d6'. 



We postpone the proof of Lemma [4.11 until the end of the present section. 
Let us denote the integral in ([HD 

7r/2 

(15) /; = J aicsm{Pi{cose))smede, 



so that evaluating the defect variance is equivalent to evaluating Ii, which 
is done in the following proposition (for even /): 

Proposition 4.2. As I — )• 00 along even numbers, we have 

(16) -^^ = + 
where 

00 

(17) Ci = j ^ (arcsin(Jo(V')) - Jo(V')) dV'- 



Moreover, the constant Ci is strictly positive. 



^T\ie asymptotic full dependence (|13|) was proven under the technical assumption 
ti,t2 7^ 0; it is natural to conjecture, though, that a slight modification of the same 
argument will work if either of ti vanishes. 
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Proof of Theorem assuming Proposition \4-^ Formula (jl4p together with 
Proposition 14.21 yields 

YariVi) = ^(1 + 0,^^(1)). 

The positivity of the constant C = 32ttCi follows directly from the positivity 
of Ci, which is claimed in Proposition 14.21 □ 



Proof of Proposition \4-^ As it was explained in Section [21 to extract the 
asymptotics of Ii, we will expand the arcsine on the RHS of (jlSp into the 
Taylor series around the origin; we will encounter only the odd moments of 
Pi{cos6), due to the arcsine being an odd function. As it was pointed out 
in Section [21 the function P/(cos0) differ from its scaling limit in that the 
integral of the former vanishes, whereas the integral of the latter diverges; all 
the other odd moments of Pi{cos6) are asymptotic to the (properly scaled) 
corresponding moments of the Bessel function, (see Lemma [5. 2p . To account 
for this discrepancy, we subtract Pi (cos 6) from the integrand and add it 
back separately. To this end we write 

arcsin(P;(cos 6*)) = (arcsin(P;(cos 0)) — P/(cos0)) + P/(c 

so that the vanishing of the integral of the latter 

7r/2 

Pi{cos 9) sin 9de = 0, 



cos I 







for even I implies 

tt/2 

(18) " / (^^csin(P,(cos6l)) - PKcos6l))sin6ld6'. 



The advantage of the latter representation (|18p over (|15p is that the only 
powers that will appear in the Taylor expansion of the arcsine on the right- 
hand side of (fTSj) are of order > 3, so that the moments of P/(cos0) are all 
identical to the corresponding moments of the scaling limit. Intuitively, this 
means that we may replace the appearances of P[{cos9) in (|18p by its the 
scaling limit. The rest of the present proof is a rigorous argument that will 
establish the latter statement. 
Let 



2k+l 

Life!. 

fc=l 



(19) arcsin(t) — t = ^ a^t 
where 

(20) - '■''y- 



are the Taylor coefficients of the arcsine. Note that all the terms in the 
expansion (jl9p are positive, and by the Stirling formula, the coefficients are 
asymptotic to 
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for some c > 0, so that, in particular, the Taylor series (jl9p is uniformly 
absolutely convergent. Therefore we may write 

7r/2 
oo L 

(22) Ii = ^ak Pi{cosef''+^sm9d6. 

k=i { 

We know from Lemma 15.21 that for every k > 1, 

Tr/2 



(23) j PiicosOf'^+^smede 



C2fc+1 



with C2fc+i given by (|33|1 : comparing this result to (j22|) . it is natural to expect 
that 

(24) = ^ + 

where 

oo 

(25) C2 = y^QfcC2fc+i- 

k=l 

We are going to formally prove ()24p immediately; however, first we eval- 
uate the constant C2 by summing up the series in ()25p . and validate that 
indeed C2 = Ci in (fT7|l . as claimed. We plug ([33]) into (f25|) to formally 
compute 

/ 00 ^ 

\2fc+l 



fc = l Q Q \fc=l 



00 



• (arcsin(Jo(V')) - Joiip)) dip = Ci, 



where to obtain the third equality, we used ()19p again. To justify the ex- 
change of the summation and integration order, we consider the finite sum- 
mation 



m ,. 



k=i 

using (j2ip and (j3ip to bound the contribution of tails, and take the limit 
m — )• 00. 

We now turn to prove (j24p . To this end we expand (jl9p into a finite degree 
Taylor polynomial while controlling the tail using Lemma [5. 71 Indeed, using 
(pT|) and (|i6]) . we easily obtain 

7r/2 7r/2 
00 00 

^ ak \Pi{cos9)\^^^^sm9d0 < ^ ak \Piicos0)f sin 0d0 

A;=m+1 Q fc=m+l q 

72 ■ 2^ 1,3/: 



1 



;2 ^3/2 /7^/2 : 

A;=m+1 ^ 
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since |-Pi(t)| < 1 for every / and t G [—1,1], so that is monotone 

decreasing with k. 

We then have for every m {m = m{l) to be chosen) 



(26) 



h = Y, ak I Piicos ef'^^^ sinOde + O (^j^) ' 

k=l n VV"^ / 



and plugging (|23p (a direct consequence of Lemma 15. 2p into ()26p finally 
yields 

(27) // = C2,m • 4 + Om ( ) + ^ 



with 

m 

C2,m = ^ akC2k+l- 
k=l 

It is clear that (j27p implies (|24p (recall the definition (|25p of C2 and note 
that as m — 7- 00, C2,m C2), which concludes the proof of the statement 
(fT6|) of the present proposition. 

It then remains to prove the positivity of the constant Ci. While its 
nonnegativity Ci > is clear since it is the leading constant for the integral 
Ii, which, up to an explicit positive constant, equals the variance of a random 
variable via (jl4p . the strict positivity is less obvious. To this end, we recall 
that Ci = C2, the latter being given by ()25p . Note that all the Taylor 
coefficients of arcsine are positive (see (j20p ). Hence the positivity of C2 
(and thus also, of Ci) follows from the second statement of Lemma 15.21 
which claims that > are nonnegative, and C3 > is explicitly given. □ 



Proof of Lemma \4-l\ The result of the present lemma is an artifact of the 
following general fact (see e.g. OdS]). Let {Xi,X2) be a 2-variate centred 
Gaussian random variable with covariance matrix 

1 r 

r 1 

\r\ < 1, and for i = 1, 2 define the random variables (recall that H{t) is the 
Heaviside function (j3|)) 



Bi = H{Xi) 



1 X, > 
-1 X- < 



Then Bi are both mean zero with covariance given by 

2 

Cov{Bi, B2) = — arcsin(r). 
vr 
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Using the definition ([3]) of the defect, we may exchange the order of taking 
the expectation and integrating to write 

(28) 

Var(PO = E[p2] =E II H{fi{x))H{fi{y))dxdy 

E[H{fi{x))H{fi{y))]dxdy = Att ^ E[H{fi{N))H{Mx))] dx, 

52x52 52 

by the isotropic property of the random field where is the northern 
pole. Note that for every x G 5^, fi{N) and fi{x) are jointly Gaussian, 
centred, with unit variance and covariance equal to 

E[fi{N) ■ fi{x)] = ri{N,x) = Piicosd), 

where {9, cj)) are the spherical coordinates of x; the latter follows from the 
definition ([2]) of the covariance function. Therefore, as it was explained 
earlier, for every x, y G 5^, 

E[H{fi{x))H{fi{y))] = -arcsin(PKcos0)). 

vr 



We then evaluate the latter integral in ()28p in the spherical coordinates as 

Var(Pi) = 8^"^ I ^ arcsin(Pi(cos 9)) sin 9 d9 = 16tt / arcsin(P;(cos0)) sin^c?^, 



which, taking into account I being even (and thus Pi{t) is also even), is the 
statement (jl4p of the present lemma. □ 

5. Moments of Legendre polynomials 

We start by recalling a basic fact on the asymptotic behaviour of Legendre 
polynomials (see for instance [H]); as usual, we shall denote by the Bessel 
function of the first kind. 

Lemma 5.1 (Hilb's asymptotics) . For any e > 0, C > we have 
(29) ^^cos9)=(^-^y\m + 1/2)9) + 6(9), 



where 

(30) d{9) < 



1/2/^3/2 

2 0<9<^' 



I 

uniformly w.r.t. I > 1, 9 £ [0,^ — e]. 

Note that we will use Lemma l5.1l onlv for the range 9 G [0, so that we 
may forget about the e altogether. We will also recall that (see again |16] ) 



(31) |Jo(V')| = O 
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Lemma 5.2. Let j > 5 or j = 3. Then 

■k/2 

(32) j Pi{cos9y sinOde = Cj^il + Ojil)), 



where the constants Cj are given by 

oo 

(33) Cj = J ijJoW^di;, 



the RHS of (|33|) being absolutely convergent for j > 5 and conditionally 
convergent for j = 3. Moreover, for every j as above, the constants cj > 
are nonnegative, and C3 > is positive, given explicitly by 

Proof. By the Hilb's asymptotics, we have 
(35) 

j Pi{cos9y smOde = j f (— ^) Mil + 1/2)9) + 6(6)] sm9de. 
^ ^ 

The contribution of the error term to ([35]) is (exploiting being bounded) 

7r/2 l/l tt/2 

(36) « J \Joiil + l/2)9)\^-^5{9)ed9 = j+ j ■ 

i/i 

Now 

i/i i/l 

\Jq{{1 + l/2)9)\^-^5{9)9d9 <^ j 9^d9 ^ ^, 

and using (|3ip . we may bound the second integral in ()36p as 

7r/2 tt/2 

^-^6(9)9d9 < , . \. / 9^/^d9 ^ 



Mil + l/2)9)r^6{9)9d9 « j^j^ J 9-y^d9 « j^—^. 

i/i i/i 

Plugging the last couple of estimates into (p6|) . and finally into ([35]) . we 

obtain for j > 3 

(37) 

I Piicos 9y Sin ed9= I (-^) Jo((/ + l/2)0)^'0d0 + o(^) . 

Therefore we are to evaluate 

■k/2 , L-k/2 

M{l^\l2)9y9dQ = -. \ .7,,,, M^)Hd^. 



sm{Q) ) ^ I ) I j2 f \^sin(V'/L) ) 
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where we denote L := / + ^ for brevity. The statement of the present lemma 
is then equivalent to 

L7r/2 .^^ ^ 

/ {Jim)' 



where cj is defined by 

The main idea is that the leading contribution is provided from the in- 
termediate range 1 <^ tp <ti eL for any e > 0; here we may replace the 
factor 

^sm{'ilj/L) J 

with 1. To make this argument precise, we write for tp G [0, ^ • L] 



s\n{ip/L) V^V ' 

so that also 

(39) f^^f + 

Therefore the integral in (j38p is 
(40) 

Liv/2 , L-k/2 



I L-k/2 
1 I' .,3, 



/ vi^o(V')P# 

V 



Note that as / — t- oo (equivalently L — t- oo), the main term on the RHS of 
converges to Cj; indeed 

Lit/2 

(41) j Jo(V)^'V'dV^c,-, 



so that it remains to bound the error term. Now 

Ltt/2 

■ ,-1 [■L-k/2 [■L-k/2 

(42) / ^^\Ui^y\d^= / + 



Jl 





and we use (|3T]) to bound the latter as 

l-k/2 



0,1)./ 



L-k/2 



SO that upon plugging the latter into ([I2]) yields 

l-k/2 

^^\jQ{^y\dij = o{i + f-^/^y 
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plugging the latter into (00]) yields 



so that (j4ip implies (|38p for j > 5, which was equivalent to the statement 
of the present lemma in this case. 

It then remains to prove the result for j = 3, for which case we have to 
work a little harder due to the conditional convergence of the integral ()33p : 
to treat this technicality we will have to exploit the oscillatory behaviour of 
the Bessel function, and not only its decay ()3ip . It is well-known that 



IT \^3/2 

SO that 

(43) /2y/^os(^-V4)^^^^ 1 



and hence the integral on the RHS (j33p is indeed convergent, by integration 
by parts. 

Now we choose a (large) parameter K ^ 1, divide the integration range 
into [0, and [K,L^]] the main contribution comes from the first term, 
whence we need to prove that the latter vanishes. Indeed, we use ()43p to 
bound 

L-k/2 



K 

where we use integration by parts with the bounded function 

T 



I{T) = j cos{tfdt 





to bound the contribution leading term in (j43p : it is easy to bound the 
contribution of the error term in (j43p using the crude estimate pip. 
On [0, K] we use ([39]) to write 



K , K 



(45) / ( Mi^f^d^P = [ Mijf^di^ + O'^"^ 



^ ^ 

the former clearly being convergent to C3. Combining (|44p with (|45p we 
obtain 

^' ^ U^fi^di; = / Ui^)Hdi; + O + > • 



K 



sin(^/L); "^^^ ^ ^ 7 ^ ^ \^ 



this implies (j38p for j = 3 (which is equivalent to the statement of the 
present lemma for j = 3) upon choosing the parameter K growing to infinity 
sufficiently slowly (e.g. K = Since j = 3 was the only case that was 
not covered earlier, this concludes proof of the moments part of Lemma [5.2[ 
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It then remains to prove the nonnegativity statement of Cj, and the ex- 
phcit expression ()34p for C3. In fact, both of those statements follow from 
the computation we performed in our previous paper [12], p. 18. We report 
the relevant results in Lemma 15.31 below. □ 

We recall an alternative characterization for moments of Legendre polyno- 
mials from |I2j. For given positive integers li,l2, h we introduce the so-called 

Clebsch-Gordan coefficients |c''^g;2o| 1 which are different from zero if and 
only if li, l2,h are such that I1 + I2 + h is even and /j + Ij < for all permu- 
tations i-iiik = 1,2,3. The Clebsch-Gordan coefficients are well-known in 
group representation theory (they intertwine alternative representations for 
5*0(3)) and in the quantum theory of angular momentum; we do not provide 
more details here, but we refer instead to [12] or to standard references such 
as [171 [18] (see also [lO]). The results provided in [12] are as follows: 

Lemma 5.3 ( [12j . Lemma A.l and above). For all even I, we have 

1 



and 

Also, for j > 5 
1 



J ' r r 



L1...L 



l...l^j_3 



Remark 5.4. The previous discussion yields the following interesting corol- 
lary: as /—)■ 00 



00 



and for j > 5 



lim 



21 



00 



L1...L 



^di). 



The following lemma establishes the lower bound for the constant C in 
Theorem 11.21 claimed in Section II. 2[ Its proof is straightforward. 

Lemma 5.5. For C as in we have 

32 



C > 



27 
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Proof. Recall that C = f Ci, where Ci = C2 is given by (p5|) . Since all the 
Qj and Cfc are nonnegative, L emma 1 5 . 5 1 follows from bounding the first term 
in the series 

32 

C > 327r X aiC3 : 



/27 

where we used ai = 1/6 and ([3l|) . □ 

Remark 5.6. As a final remark, we note that the moments of Pi{cos9) 
are themselves of some physical interest. In particular, while analyzing 
the relationship between asymptotic Gaussianity and ergodicity of isotropic 
spherical random fields, [TlJ considered the random field 



T{x) ■.= J2Tiix) = J2ciYio{9 



1=1 1=1 



where Yio{9, = \J -^Pi{cos9) denotes standard spherical harmonics (for 

m = 0) and g € 5*0(3) is a uniformly distributed random rotation in M^. 
It is readily seen that the resulting field is isotropic, and each of its Fourier 
components at frequency / is marginally distributed as (for any x ^ S"^) 



Yi = ci^^^^Pi{t),t^U[0, 1]. 

Up to normalization constants, then, the moments of the marginal law are 
exactly those of the Legendre polynomials we established earlier in the paper. 
More precisely, if we focus (as in [11]) on Ti := Ti{x)/ y^Var{Ti), we obtain 
immediately, for any x G 5^ 

E\fi{x)} =E{^{2l + l)Pi{t)} =0, e{t)(x)}' = 1, 

and 

E {fiix)}' = S {^/(27Tl)PK*)}' = o (^^) , 

whereas J2] 

ii;{V(2/ + i)PKi)}'-iog/ , 

and hence all moments of order q > 4 diverge. 

Note that in Lemma [5. 2 1 we worked relatively hard to establish the precise 
asymptotics for the moments of Legendre polynomials. A much cruder ver- 
sion of the same argument gives a uniform upper bound for the 5th moment 
of the absolute value of the Legendre polynomials: 

Lemma 5.7. We have the following uniform upper bound for the 5th mo- 
ment of the absolute value of Legendre polynomials 

Tv/2 

(46) / \Pi{cose)fsmede = o(^ 



where the constant involved in the ^O' -notation is universal. 
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Proof. The proof follows along the same lines as the beginning of proof 
of Lemma 15.21 for A; > 5, except that we use the crude upper bound (j3T]) 
whenever we reach (j37p . and the trivial inequality 



< 1, for e G 



sin( 



0, 



vr 



□ 
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